Josephson effect due to the long-range odd-frequency triplet superconductivity 

in SFS junctions with Neel domain walls 
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We consider a SFS Josephson junction made of two superconductors S and a multidomain fer- 
romagnet F with an in-plane magnetization. We assume that the neighboring domains of the 
ferromagnet are separated by Neel domain walls. An odd-frequency triplet long-range component 
of superconducting correlations arises in the domain walls and spreads into the domains over a long 
distance of the order £t = \/D/2tyT, where D is the diffusion coefficient (dirty limit is implied). 
We calculate the contribution of this component to the Josephson current in the situation when 
conventional short-range components exponentially decay over the thickness of the F layer and can 
be neglected. In the limit when the thickness of the F layer is much smaller than the penetration 
length of the long-range component, we find that the junction is in the tt state. We also analyze a 
correction to the density of states due to the long-range triplet component. 

PACS numbers: 74.50.+r, 74.45.+C, 74.78.Fk, 75.70.Kw 



I. INTRODUCTION 

The past decade was marked by a rapid growth of in- 
terest to the study of hybrid superconductor-ferromagnet 
(SF) structures (see, for example, reviews [31211) . New 
physical phenomena arising in these systems originate 
from a nontrivial interplay of competing orders in su- 
perconductors and ferromagnets. On one side, electron- 
electron interactions lead in superconductors to forma- 
tion of Cooper pairs consisting of two electrons with op- 
posite spins. On the other side, the exchange interaction 
in ferromagnets tends to align the electron spins parallel 
to each other. In SF structures these two types of inter- 
actions are spatially separated and can coexist despite 
much greater value of the exchange energy h in compar- 
ison with the superconducting gap A. 

Due to the proximity effect^ the superconducting cor- 
relations penetrate into the ferromagnet in SF structures. 
The condensate wave function / penetrates into the ferro- 
magnet with a uniform magnetization M over a distance 
of the order of the "exchange length" £/, = y/D/h. 1,2 
The condensate wave function decays in F in a nonmono- 
tonic way as f(x) ~ exp(— x/^) cos(x/£h): it oscillates 
in space and decreases exponentially. This nonmonotonic 
behavior of /(x) leads, in particular, to a it state in SFS 
Josephson j unctiona ^ 2 »L S i Z i 2i2il2i characterized by a neg- 
ative critical current I c in the Josephson current-phase 
relation I (if) = I c siny. 

If the magnetization in the ferromagnet is nonuniform 
a new phenomenon becomes possible: a triplet compo- 
nent of the condensate wave function / (generally speak- 
ing, the condensate wave function is a matrix in the 
particle-hole and spin space) arises in the SF system, 3 
This triplet component is an odd function of the Mat- 
subara frequency w (while the conventional BCS singlet 
component of / is an even function of lo) and spreads 
in the ferromagnet over a long distance of the order of 



£t = \/D/2ttT. The existence of this long-range triplet 
odd-frequency component in SF structures with an inho- 
mogeneous magnetization was predicted in Ref. [ll| and 
further discussed in Ref. 12. Unlike the triplet component 
in superfiuid 3 He and in Sr 2 Ru04, this odd-frequency 
triplet component corresponds to s-wave pairing and, 
hence, is symmetric in the momentum space. Therefore 
it is not destroyed by scattering on nonmagnetic impuri- 
ties and survives in the dirty limit. We call this compo- 
nent the long-range triplet component (LRTC). Histori- 
cally, the odd-frequency triplet pairing was conjectured 
in 1974 by Berezinskiii^ as a possible mechanism for su- 
perfluidity in 3 He but this conjecture was not confirmed 
experimentally. 

There is a significant amount of experimental data 
that may be interpreted as manifestation of the LRTC 
in SF systems^^^iii^^^iSi^ Of particular impor- 
tance are the experiments on SFS systems in which a 
long-range phase coherence of the condensate wave func- 
tions was observed in ferromagnets with a length con- 
siderably exceeding the penetration length of the singlet 
component £,h' 1S]21]22 In principle, this long-range phase 
coherence can be due to the LRTC but, for unambigu- 
ous identification of the LRTC, further experimental and 
theoretical studies are very important. In particular, the 
theory of the LRTC was developed for specific types of 
the inhomogeneous magnetization in ferromagnets, which 
do not exhaust all possible types of the magnetic struc- 
tures in real samples. 

In Refs. H0 and |H the LRTC was studied in SF 
systems with a Bloch-type magnetic structure [the mag- 
netization vector M(x) lies in the y-z plane parallel to 
the SF interface and rotates with increasing x] . The am- 
plitude and the penetration length of the LRTC induced 
in the ferromagnet have been calculated in Refs. [ll| and 
HH - Under certain conditions the LRTC penetrates the 
ferromagnet over a long distance of order As shown in 
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Ref. HH, where a SFS structure with a conical ferromag- 
net was studied, the LRTC may decay in a nonmonotonic 
way provided the cone angle exceeds a certain value. In 
Refs. [24] and [2^ a multilayered SF structure was investi- 
gated in which the magnetization vector M(cc) has fixed 
but different orientations in different ferromagnetic lay- 
ers. This structure is also similar to a Bloch-type domain 
structure. 

At the same time, it is well known that the domain 
structure in a ferromagnet can be not only of the Bloch 
type but also of the Neel type^ A SF system with a 
Neel- type magnetic spiral structure [the M(y) vector lies 
in the y-z plane parallel to the SF interface and rotates 
with increasing y] was studied theoretically in Ref. [27l . 
This magnetic structure may be regarded as an infinite 
Neel wall. It was shown that in this case the LRTC does 
not arise in the system. However, this statement is valid 
only for the case of a uniformly rotating M(y) vector. In 
a more realistic situation of a ferromagnet with magnetic 
domains separated by Neel walls, the LRTC arises at the 
domain walls and decays inside the domains over a long 
distanced Another source of the triplet component in 
SF structures was considered in Ref. |29|, where the SF 
interface was assumed to be spin-active. The Josephson 
effect in a S/HM/S junction (HM is a ferromagnetic half- 
metal) was studied in this paper and it was shown that 
the critical Josephson current has a maximum at low, but 
nonzero, temperature. 

In the present paper we consider a SFS Josephson junc- 
tion with the Neel domain structure in the ferromagnetic 
layer and calculate the Josephson critical current I c in 
this system. Up to now there have been no theoretical 
investigations of this problem. In Refs. H3, H3, HI, and 
I23I the critical current I c was calculated for a magnetic 
structure similar to the Bloch type. A SFS Josephson 
junction with a rotating M(x) and with the thickness 
of the F film d of order was studied in Ref. [3(3- It 
was shown that the rotation of the magnetization vector 
M(a;) leads to the appearance of the LRTC and to a sup- 
pression of the 7r state in the Josephson junction. The 
situation in a multilayered SF structure with noncollincar 
magnetizations in the F layers is more complicated. The 
sign of I c depends on the chirality, that is, on whether 
the M vector rotates in the same direction in the whole 
system or it oscillates with respect to the z axisi 24 i 25 In 
a SFS junction with a conical ferromagnet of thickness d 
much greater than £/, , the critical current I c is due to the 
LRTC and the sign of I c depends on the cone angled 

Note that in ferromagnets with a multidomain struc- 
ture the LRTC does not arise if the thickness of the do- 
main walls is very small. In this case the magnetization 
vectors M are collinear and only singlet and short-range 
triplet components are induced due to the proximity ef- 
fect. The critical current in a S FS j unction with such a 
structure was calculated in Ref. [3l|. The domain struc- 
ture also results in a suppression of the ir state due to an 
effective averaging of the exchange field. Moreover, if the 
exchange field h changes its sign over a scale shorter than 



max(^/i, I) (where I is the mean free path), long-range ef- 
fects arise even in the absence of the LRTC (the exchange 
field is effectively averaged out). This case is realized 
provided an antiferromagnet (AF) is used in SFS junc- 
tions instead of a ferromagnet. The Josephson current in 
S/AF/S junctions was studied theoretically in Refs. 
and [HI and experimentally (see Ref. [34] and references 
therein). 

We consider a domain structure in a thin F film, where 
domains with antiparallel in-plane magnetizations are 
separated by the Neel walls (while the magnetization 
does not change across the thin F film). This domain 
structure is realized in real ferromagnetic filmsi 2 ^ The 
y-z plane is chosen to be parallel to the SF interfaces 
(see Fig. [T]). We show that the LRTC arises at the Neel 
domain walls and decays exponentially away from the 
domain walls and the SF interfaces over a long distance 
£t- We calculate the Josephson current due to the LRTC 
and find that its sign corresponds to the 7r junction. The 
mechanism of the n junction in our case is related to 7r/2 
phase shifts at the SF interfaces, while the LRTC does 
not oscillate inside the F layer (in contrast to the short- 
range component). We also study modifications of the 
density of states in the F film due to the LRTC. 

The paper is organized as follows. In Sec. [HI we formu- 
late the Usadel equations and the corresponding bound- 
ary conditions, investigate the main features of the long- 
range triplet superconducting component that appears 
due to the presence of Neel domain walls, and find the 
Josephson current due to this component. The analysis is 
made for the simplified model with only one half-infinite 
region with rotating magnetization. In Sec. IIII1 we con- 
sider the realistic case of the multidomain F layer using 
the results of Sec. [TIJ In Sec. IIVI we study the correc- 
tion to the density of states due to the triplet component. 
This section is an addition to our previous paper Ref. |28| . 
Finally, we discuss limitations of our model in Sec. [V] and 
present our conclusions in Sec. IVII 



II. BASIC EQUATIONS. JOSEPHSON EFFECT 
IN THE CASE OF A HALF-INFINITE DOMAIN 

We consider a ferromagnet (0 < x < 2d) sandwiched 
between two bulk superconductors. We assume that the 
in-plane exchange field h(y) in the F layer proportional 
to the magnetization M, depends on y as follows: h(y) = 
h(0, 0, 1) at y < and h(y) — h(0, sina(y), cosa(y)) with 
a {y) = Qy at y > 0. This means that the magnetization 
vector M is oriented along the z axis at y < and rotates 
in the y-z plane at y > 0. The region with the rotating 
magnetization models the Neel domain wall. This struc- 
ture is shown in Fig. [lja) and contains only one half- 
infinite domain and one half-infinite region with rotating 
magnetization. Then we shall use the obtained results 
to describe a realistic multidomain structure depicted in 

Fig. him. 

In order to calculate the Josephson current, we need to 
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FIG. 1: SFS junctions considered in the paper, (a) The do- 
main (y < 0) and the region with rotating magnetization 
(y > 0) in the F layer are half-infinite, (b) Multidomain F 
layer. Depending on the relative orientation of rotating mag- 
netizations in the neighboring domain walls, we distinguish 
the cases of positive and negative chirality (Q has the same or 
opposite sign in the neighboring domain walls, respectively). 



find the condensate Green functions in the ferromagnet 
induced due to the proximity effect. We consider the 
dirty limit, which means, in particular, that hr <C 1, 
where r is the momentum relaxation time due to elastic 
scattering. 

In the dirty limit, the system is described by the Us- 
adel equation for the matrix Green function g which is a 
4x4 matrix in the Gor'kov-Nambu and spin spaces. The 
Usadel equation for the case of an inhomogeneous mag- 
netization was written in Refs. [35| and y. However, we 
redefine the Green function used in Refs. HH and[H (and 
also in our previous paper Ref. l28h introducing a new 
matrix function g ncw — VgV^ with the transformation 
matrix 



V = exp 



('!< 



i~3 - T )a 3 



(1) 



where fj and <ii are the Pauli matrices in the Gor'kov- 
Nambu and spin spaces, respectively. This transforma- 
tion was proposed in Ref. l36l . and below we shall use 
the new Green function g ncw omitting the subscript for 



brevity. The convenience of the new definition is that 
the Usadel equation for the new Green function possesses 
the explicit symmetry with respect to rotations of the ex- 
change field h: 

DV(gVg)-u[T 3 Z ,g]-i[T 3 (hcr),g}- [A,g] = 0, (2) 
A = (n Re A - f 2 Im A) a , 

where u> is the Matsubara frequency, h = hn, n = 
(0, sina(y), cos a(y)) and a(y) = at y < while 
a (y) — Qy at V > 0- We assume that the diffusion 
coefficients D for electrons with spins up and down are 
equal to each other (this is correct provided the exchange 
energy h is much less than the Fermi energy sp). 

The Usadel equation ^ is written in the general form 
taking into account both the superconductivity and mag- 
netism. In the SFS junction, the pair potential A is 
nonzero only in the S layers, while the exchange field 
h is nonzero only in the F layer. There is no attrac- 
tive interaction between electrons in the F layer, hence 
A = 0. However, the condensate (Gor'kov) functions are 
finite in the F region due to the boundary conditions at 
the SF interface. 

In the general case, the Green function has the follow- 
ing components: 

9 = ^3 (.go^o + g<r) + fi (/oo-o + fo") + f 2 (fo&o + f<r) . 

(3) 

In the bulk of a normal metal, only the go component is 
present. The superconducting correlations are described 
by the scalar anomalous components /o and /o, while 
the vector components (g, f , f) are due to the ferromag- 
netism. 

The superconducting correlations described by the / 
components (nondiagonal in the Nambu-Gor'kov space) 
are assumed to be weak due to a finite interface trans- 
parency. In the considered case of weak proximity effect 
(|/| <C 1), the Green function acquires the form 



g = f 3 <T sgnw + F, 
where the anomalous part can be written as 



(4) 



(5) 



F = T 1 f + T 2 f 

with matrices in the spin space 

/ = f ao + {&, (6) 
f = f a + {&. (7) 

Equation ([2]) can be linearized and brought to the form 

V 2 F-2klF-ik 2 h sgnu{T (™),F} = 0, (8) 

where u = nT(2n + 1), fc£ = \u\/D, k\ = h/D, and the 
braces denote the anticommutator. 

The anomalous part of the Green function in the bulk 
of the superconductor with the superconducting phase <p 
is Fs{<p) — (ti cos tp — T2 sin <p)crafs, where 



fs 



Vuj 2 + A 2 



(9) 
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We intend to find the Josephson current at the phase dif- 
ference (p between the two superconducting banks. As- 
suming the phases of the left and right superconductors 
to be — ip/2 and ip/2, respectively, we write the boundary 
conditions for F at the SF interfaces: 



dF 

dx 



lb 



(10) 



where the two signs correspond to the left [x = 0) and 
right (x — 2d) SF interfaces, respectively. Here -fb = Rb&, 
while Rf, is the interface resistance per unit area and a 
is the conductivity of the ferromagnet. This boundary 
condition follows from the Kupriyanov-Lukichev ones^l 
provided two assumptions are satisfied: (1) the proximity 
effect is weak (i.e., Jb/^h 3> 1) and (2) the bulk solution 
Fg in the superconductors is unperturbed and valid up to 
the interface (i.e., jb/^g S> cr/<Jg, where £g — yjDg/A 
and ag are the coherence length and the conductivity in 
the S banks). 

The technical problem with Eq. ([5]) is that this is a 
two-dimensional partial-derivative differential equation. 
However, we can employ a trick similar to the one pro- 
posed in Ref. |28|, which allows us to make the Fourier 
transformation over x, reducing the problem to only one 
dimension y. 

The equations for the / and / functions split and for 
the / function we obtain 

V 2 f-2klf-ik 2 h sgnu J [{na) 1 f] = 0, (11) 

(12) 



dx 



T 



fg cos f 



76 



<Tq. 



The / function obeys the same equation, although the 
boundary conditions are different: 



V 2 f-2k 2 J-ik 2 h sgnu J [{n*)J} 



df 
dx 



x=0,2d 



fs sin I 

76 



(T - 



(13) 
(14) 



The functions / and / are defined for < x ^ 2d. 
The function / is even with respect to the center of the 
F layer, while / is odd. We can continue the functions to 

the whole x axis: for / we do it periodically and for / — 
antiperiodically, obtaining continuous functions in both 
cases. Due to the boundary conditions both the functions 
have cusps at x = 2dN with integer N . The boundary 
conditions (|10j) producing the cusps can be incorporated 
into Eq. © itself as 5-functional terms (this situation 
is similar to the standard quantum-mechanical problem 
with ^-functional potential producing a cusp of the wave 



function): 



V a /-2A2/-t*2 S gaa;{(nd-),/} 



-O-Q 



C\ P (f CO 

2fsco Sl ^ 6 (x-2dN), (15) 



76 



N- 



V 2 f-2k 2 J-ik 2 sgnu J {(n&)j} 



-0"o 



/sSUlj J2 (-lfS(x-2dN). (16) 



76 



N- 



Now, instead of solving Eq. ((5|) at ^ x ^ 2d with the 
boundary conditions (fTT)]) . we have to solve Eqs. (| 1 5[) and 
([H]) at all x. 

The Fourier transformation 

f(k,y)= dxexp(-ikx)f(x,y) (17) 



in Eqs. ([15]) and (|16|) should be performed over "bosonic" 
wave vectors k n and "fermionic" wave vectors k n for the 
periodic function / and antiperiodic function /, respec- 
tively: 



7T - 7T 

k n = — 2n, k„ — — (2n + 1). 
2d 1 " 2d ' 



(18) 



We obtain 



2 is cos ■ 



-^1 - (kl + 2k 2 )f - ikl sgrnu {(no-), /} 

'"(19) 

g - (kl + 2kl)f-ikl sgnc {(n«x), /} = 2 

"(20) 



76 



-co, 



a a. 



The two equations are similar and we may consider 
only one of them, say Eq. (I19|) for the / function. Then 

the result for the / function can immediately be obtained 
by substituting k n <— > k n and cos ^ ^ sin ^ . 

At y > the function a(y) is y-dependent, while at 
y < we have a = 0. In the region of positive y one can 
exclude the y-dependence from Eq. (fT^|) with the help of 
a rotation 

f = Uf lt U + , (21) 
where U = exp(i(Tia(j/)/2) . As a result, we get (y > 0) 



' J2fu ( kl + % + 2k 2 \ h+^-hfoi+iQi*!, ^] 



dy 2 



ikl sgnw{a 3 , /„} 



2/g cos | 

76 



dy 
a (22) 
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in terms of the new function f u (k, y) (the square brackets 
denote the commutator). The same equation is valid for 
y < if we put Q = 0: 



2 fs cos £ 



76 



CT . 



<9y 2 

(23) 

The original functions / and df jdy are continuous at 
y = 0. Therefore, the rotated functions obey the follow- 
ing boundary conditions at y = 0: 



A.(-o) - /„(+o), 

S/f(-0) _ 9/ u (+0) 



<9y 



dy 



+ i 



Q 



(24) 
(25) 



Thus, we have to solve the linear matrix Eqs. (|22fl (y > 
0) and (f2"3"| (y < 0) of the second order with the boundary 
conditions ([24]) and (|25[) at y = 0. We can represent the 
solution in the form 



(26) 



where 9 is the Heaviside step function and the constants 
fi(Q) and F(0) are the homogeneous solutions of Eqs. 
and (|23p at y = ±oo. The matrices T have the form 



where 



2/ s cos f (k 2 n + Q 2 + 2fc 2 



^a(Q) = - 



4^ / g cos f k\ sgn 
7fW) 



(27) 

(28) 
(29) 



and 



V(Q) = (kl + Q 2 + 2kl) {hi + 2kl) + Akt (30) 

The correction 6f u (k,y) obeys the same Eqs. (|22p and 
without the right-hand side. It has the form 



Sf u = <Wo + 0-3/3 + 0-2/2- 



(31) 



The first term in Eq. (|3Tj) is the singlet component. The 
second term is the triplet component with zero z projec- 
tion of the Cooper-pair spin. This component arises even 
in the case of a homogenous magnetization of the ferro- 
magnet and decays in the F film at the short distance 
The last term in Eq. (|31[) is the triplet component with 
the spin projection ±1. It arises in the case of an inhomo- 
geneous magnetization and decays over a long distance 
of the order £p. The functions fi(k,y) in Eq. (|3"Tj) can be 
represented as a sum of eigenf unctions of Eqs. (f2"2")l and 



My) =J2 A il e M~ K l(Q)y) aty>0, (32) 
I 

fi(y)=J2 B » e M K i(°)y) aty<0. (33) 



The inverse decay lengths Ki(Q) are the eigenvalues of 
Eqs. (H2J) and d23J) (without the right-hand side). The 
equation for Ki(Q) has the form (I — 1, 2, 3) 

(nf - kl -Q 2 - 2fc 2 ) 2 + 4 [Qmf] (« 2 - kl ~ 2k 2 J) 

+ 4fc,t( K 2 -fc 2 -Q 2 -2fc 2 ) =0. (34) 

We assume that the exchange length is the shortest 
length in the problem: 



k h ~^> k n ,Q ,k u 



(35) 



Then, the eigenvalues ki consist of two "short-range" val- 
ues 



«± ~ (1 T i sgn uj) kh 
and one "long-range" value 

K L [Q) « y/kl + Q 2 + 2kl 



(36) 



(37) 



At y < we have the same k; with Q = 0. 

Calculating the corresponding eigenvectors and match- 
ing solutions (|32|) and (I33| with the help of the boundary 
conditions (|24|l and (f25|) . we find the coefficients Ay and 
£?;z . This simple but cumbersome calculation is similar to 
the one presented in Ref. 28. In the considered limit of a 
small exchange length [see Eq. the coefficients A 2 l 

and B2L that describe the LRTC are the largest ones, 



.4 



2L 



B- 



2L 



QF3 



(38) 



where for brevity we have denoted 

kq = kl{Q), k q = k l (0). (39) 
In the limit (|3"5|) the function Tj,(Q) has a simple form 



^3 



^/s cos I sgn a; 
lbk 2 h 



(40) 



Therefore, the magnitude of the LRTC at the interface 
between the domain and the domain wall (y = 0) is equal 
to 



f L (k ni 0) = f 2L (k n ,0) 



i/scosfsgnw Q 



%ki 



KQ + Kq 



(41) 

while the decay along the y axis is determined by kq 
or kq at y < and y > 0, respectively [see Eq. ([52])]. 
Below we consider the situations when only the LRTC is 
essential, and denote the corresponding contribution /jl 
to the Green function by /l for brevity. 

The real-space function is determined by the inverse 
Fourier transform 



2d 



e l ^ x f L (k n ,y). 



(42) 



The /t function is obtained after substituting fc n 1— > fc„ 
and cos $ 1— > sin ^ . 
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A. Josephson current 

The supercurrent density is determined by the anoma- 
lous part of the Green function: 



(43) 



where a is the conductivity in the normal state (we shall 
calculate this expression inside the F layer, hence a is 
the conductivity of the F material) . Since we choose the 
phases of the left and right superconductor as then 
the / components are even with respect to the center of 
the interlayer, while the / components are odd. Calcu- 
lating the supercurrent in the center of the F layer, we 
obtain 



Tr (f 3 a„FVf) = 4z (/ V/o + fVf) 



(44) 



Only the x component of the current survives due to the 
odd behavior of the / functions. 

However, calculating the current not in the center of 
the F layer we also obtain the y component. In the region 
of the F layer, where only the long-range component is 
present, we have F = tio^/l + t 2 (T2 II an d finally 



(45) 



In the region with the magnetization rotating as a func- 
tion of y, the final expression remains the same with fi 
and Jl being the second components of /„, see Eqs. ([2"Tj). 
(HI, and (j3T|). 

Inside the half-infinite domain, Eq. (|42p immediately 
yields 



fh{x,y) =~y^e 



ifs cos(ip/2)Q sgnw 
2d^bkl(K Q + k ) 



(46) 



M*,y) = - 2Je ^yffi 2 > ^ e-ov, (47 ) 
^ 2d<y b kf l (K Q + K Q ) 



In the limit of thin F layer, i.e., d <C 1/Q.l/k^ (at 
the same time, we assume 1/kh <C d, which allows us to 
neglect the short-range components), we obtain 



fh{x,y) = - 



ifs cos(ip/2) sgnw 



2d lb k 2 



VQ 2 + 2kl + V2K 



fe„>0 



^ Q- cos k n x e knV j , 
(48) 



(49) 

We calculate the Josephson current according to Eq. 
(I45p . In the limit of thin F layer, the main contribution 



is given by the second term in the brackets /lV/l, more- 
over, it is sufficient to keep only the k n = component 
in Jl [the first term in the brackets in Eq. (|48|) ]. This 
is because the components containing k n in the denom- 
inator in Eq. (|48[) are much smaller. Summing over k n 
finally yields 



i{x,y) = - simp 



aQ 2 irT 
\Qe{d lb k 2 ) 2 



E 



/I 



■ exp 



VQ 2 + 2kl + V2K 
e x cosh (f§) sin (ff ) + e y sinh (%) cos (ff ) 



sinh 2 (|f)+sin 2 (| 



(50) 



2,/ , 



where e x and e y are the unit vectors in the x and y 
directions, respectively. In the region with Q ^ 0, we 
obtain Eq. ^ with ^JQ 2 + 2k 2 instead of y^k 2 in the 
argument of the exponential. The main qualitative result 
of this formula is that the x component of the current 
density has a form j x — j cx simp with negative j cx , i.e., 
the 7r state of the junction is realized. 

Expression (|5D|) for the supercurrent is valid in the re- 
gion where only the LRTC component is essential, while 
the short-range components are exponentially small. 
Therefore, close to the SF interfaces and to the boundary 
between the domain and the domain wall (at distances 
of the order of £/,) the expression is not applicable. 

Note that Vj = and V x j = within our accuracy 
(we must neglect k^ and Q in comparison with 1/d). 
The total y current is absent, while the total x current is 
equally shared between the regions with constant and ro- 
tating magnetization and corresponds to the tt junction: 



2d 

j y (x,y)dx = 0, 

roo 

j x (x,y)dy = / j x (x,y)dy 

-oo JO 



16ed 7 2 fc, 



VQ 2 + 2fc 2 + V2k u 



(51) 



(52) 



The results for the x component of the critical cur- 
rent density j C x{x,y) are illustrated in Fig. [21 where we 
plot \j cx \ as a function of y at several values of x. The 
current density is maximal at the boundary between the 
domain and the domain wall (y = 0), and the maximal 
current across the junction is carried along this line. At 
the same time, the relative height of this peak depends 
on x. The distribution is smoothest at the center of the 
F layer (i.e., at x — d), while the peak at the domain- 
domain-wall boundary becomes more pronounced closer 
to the SF interfaces. 

In order to estimate the absolute value of the critical 
Josephson current following from Eq. (|52p , we assume 
that the junction area is 50 x 50 jitm 2 , a ~ (50 cm) -1 , 
Q ~ (50 nm)-\ d ~ 7 6 ~ 56m h ~ 500 K, and D ~ 
10 cm 2 /s. Then the critical current at low temperatures 
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d/2 

d (middle of 
the F layer) 




FIG. 2: SFS junction with a half-infinite domain: absolute 
value of the critical current density j cx in the x direction 
as a function of the y coordinate at low temperatures. The 
current is negative, i.e., the tv state of the junction is realized. 
The normalization constant is jo = oDQ 3, /8e(d'ytkf l ) 2 . Other 
parameters: DQ 2 = A and Qd = 0.1. 



is of order 3 /zA which is well within the experimentally 
measurable range. 



III. JOSEPHSON EFFECT IN A 
MULTIDOMAIN SFS JUNCTION 

In this section we study the LRTC in a SF structure 
with a multidomain ferromagnetic layer [Fig.[IJb)]. One 
can distinguish between two possibilities: (a) positive chi- 
rality, when the magnetization vector M(y) in all the 
domain walls rotates in the same direction (e.g., clock- 
wise), and (b) negative chirality, when the vector M(y) 
in neighboring domain walls rotates in the opposite di- 
rections [e.g., clockwise in the 2nth domain walls and 
counterclockwise in the (2n + l)th domain walls]. We 
are interested in the LRTC assuming that the exchange 
length £/j is much smaller than the coherence length £y. 
At distances x essentially exceeding the length £/, only 
the LRTC survives in the F layer. 

We assume that the width of the domains with Q = 
is 2arj and the width of the domain walls (Q ^ 0) is 2oq. 
The origin (y = 0) is located in the middle of a domain 
with the constant magnetization. At x 3> £h only the 
long-range components of the condensate function sur- 
vive in the ferromagnet. The largest long-range compo- 
nent is the LRTC. At the boundary between a domain 
and a domain wall the solution must satisfy boundary 
conditions (J2U) and (|2"5), 



A. Positive chirality 



is also odd — this general symmetry can be demonstrated 
in Eq. ([T5J). Hence the solution for the LRTC is 

f L (k n ,y) = Asmh(K y), -a a <y<a , (53) 

hikniV) = Bsmh(K Q (y - a - a Q j), a < y < a + 2a Q . 

(54) 

Matching these solutions and their derivatives at y = ao> 
we find 



D 



-A 



sinh(KoOo) 
sinh(nQaQ) 



COSh(KQOQ) (^KQ + KQ 



tanh(KQaQ) 
tanh(Koao) 



The amplitude of the LRTC at y = ag is 



fh{k n ,a ) 



kq coth(KQOQ) + k coth(ft ao) 



(55) 



(56) 



We see that /i(fc ra ,ao) turns to zero both at clq — > 
and ao — > 0. These limits mean that the widths of the do- 
main walls and domains are assumed to be small in com- 
parison with £t but larger than The case clq = im- 
plies that we have a domain structure with the collinear 
orientation of magnetizations. The case ao = corre- 
sponds to a SF structure with continuously rotating mag- 
netization (the case studied in Ref . l27t) . In both the cases, 
the LRTC does not arise. 

The spatial dependence of the LRTC in the domain 
(\y\ < ao), corresponding to Eq. ([56]) . is given by the 
inverse Fourier transformation 



h(x,y) 



1 

2d 



E 



fL(k„,a 



sinh(K y) 
sinh(/to&o) 



(57) 



Interestingly, the function /l(x, y) turns to zero in the 
center of the domain (y = 0). This means that the 
Josephson current due to the LRTC also turns to zero 
in the domain center. 

Formula (f57|) can be drastically simplified in the limit 
when the F film is thin for the long-range component but 
thick for the short-range one (i.e., kh ^ 1/d ^> Q : k u ). 
In this case, the main contribution is given by the n = 
harmonic (with k n = 0) , since otherwise kq and kq in the 
denominator of Eq. I|56p become very large. Employing 
also Eq. (|40[) we therefore obtain 

ifs cos(y/2) sgn uj Q sinh (V^k^y) 
JL\y — r- — p2 ~ ■!//?:, 7' y* 3 ") 

where we have denoted 



/C+ sinh (\/2fcw<xo) 



sjQ 1 + 2kl coth(v'Q 2 + 2kla Q ) 
+ V2k u coth(V2fc w a ). 



(59) 



Consider first the case of positive chirality. The angle 
a(y) is then an odd function of y, which means that f2(y) 



The x dependence has vanished since the F layer is thin 
and the even function /l(x) is nearly constant. 
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Now we can easily write down the result for fa, em- 
ploying the rule formulated in Sec. [TTJ in Eq. (|57|) we 
should substitute k n i— > fc„ and cos ^ i— > sin § . After 
that we make the final step, assuming the limit of thin 
F layer. This step is different from the case of the Jl 
function, because there is no mode with k n = 0, hence 
we must retain all the modes in the sum over n: 



E 



^/s sin(y/2) sgncj 
2d lb k 2 
2Q cos k n x 



sinh (k n y) 



k n > 



' k n [coth (k n a Q ) + coth (k n a )] sinh (fc„a ) 



(60) 



Finally, we find the supercurrent: 



<rQ 2 nT x ^ /J sinh(V2^2/) 



4e(d 7& fc 2 ) 2 /c+sinh(y2fc w a ) 
E e lS ,„ ft „) s ,„ h ( t ,„)-e,co 3 ft.,)co S M^) 



k n >0 



[coth(fe„aQ) + coth(fc„ao)] sinh(/c„ao) 



(61) 



If we want to calculate the current in the region with 
Q =/= 0, then we obtain the same expression, but with 
y counted from the center of the domain wall, with 
\JQ 2 + 2k 2 instead of y / 2fcj in the arguments of two 
sinh's, and with aQ instead of ao in the arguments of two 
sinh's. For example, at a < y < a + 2oq, we have 



. . oQVT J|sinh(VQ 2 + 2fc2y' 



4e(d 7b fc2)2 ^ /c+ s i nh(v /Q2 + 2fc 2 aQ) 



fe„>0 



CO 



th(fc n aQ) + coth(fc n a )] sinh(fc„aQ) 



(62) 



where y' = y - a Q - a Q . 



From the boundary conditions (|24[) and (|25|) we find the 
coefficients C and D, and finally obtain 



fL(k n , ao) = 



kq coth(KQOQ) + kq tanh(K ao) ' 



(65) 



In this case the LRTC disappears only in the limit aQ — > 
because in this limit one again has a domain structure 
with the collinear orientation of magnetizations and very 
narrow domain walls. 

Further analysis is similar to that for the previous case 
of positive chirality. Inside of the domain (\y\ < ao) we 
obtain 



h(x,y) = ^^2e ik ^ x f L {k ni ao\ 



cosh(KQ^) 
cosh(Koao) 



(66) 



In the limit k h > 1/d >Q, k u , Eq. ([66]) yields 
if s cos{ip/2)sgnuj Q cosh (V^k^y) 

jL(y) = 77i — n t- n~7K, 7> ( b7 ) 



h(x,y) 

*E 



2dj b kf l K,- cosh (V2fc w a ) 
if s sm(tp/2) sgnw 



2d lb k 2 
2Q cos k n x 



cosh (k n y) 



k„>0 



k n [coth {kn a o) + tanh (fe n ao)l cosh {k n ao) ' 

(68) 

where we have denoted 



/C_ = v/Q 2 + 2kl coth(v/Q 2 + 2k 2 u a Q ) 

+V2k UJ tanh(V2k UJ ao). (69) 
Finally, we find the supercurrent: 

aQ 2 nT ^ f 2 s cosh(V2fc w y) 



jo(a:,y) = - sin yj- 



E 



E 

fc„>0 



4e{d-y b k 2 ) 2 ^ K_ cosh( N /2fc (i) a ) 
sin(fc„x) cosh(fc ra ?/) — cos(fc„a;) sinh(fc„y) 



yfch(fc n aq) + tanh(fc„ao)] cosh(fc„ao) 



(70) 



Similarly, in the region with Q ^ (at ao < ?/ < &o 
2aq>), we obtain 



B. Negative chirality 

Consider now the case of negative chirality, when the 
M(y) vector rotates in the opposite directions in neigh- 
boring domain walls. In this case the spatial dependence 
of /^(y) in the domain walls remains the same as before, 
i.e., this function is odd with respect to the center of 
a domain wall. However, the spatial dependence of the 
LRTC in the domains changes drastically: it becomes 
an even function with respect to the center of a domain. 
Therefore this dependence is 

fhikn.y) = Ccosh(K y), -a <y<a a , (63) 

fa{k n ,y) = Dsmh(K Q (y - a - a Q )), a < y < a + 2a Q . 

(64) 



. aQ^rT x . f 2 smh(^/Q 2 + 2k 2 ,y>) 



E 

k n >0 



4e{d lb k 2 ) 2 ^ JC- smh(^/Q 2 + 2k 2 a Q ) 
e x sin(fc„x) cosh(k n y') — e y cos(k n x) sinh(fc„y') 



[coth(fc„aQ) + tanh(fc n a )] sinh(fc„ag) 



(71) 



where y' — y - aQ - a . 



C. Discussion of the results 

The expressions obtained for the supercurrent are valid 
in the region where only the LRTC component is essen- 
tial, while the short-range components are exponentially 
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a oQ 

FIG. 3: Multidomain SFS junction: the absolute values of 
the critical currents J c + (positive chirality) and J c _ (nega- 
tive chirality) per period of the domain structure as functions 
of the domain half-width ao at low temperatures. The cur- 
rents are negative, i.e., the ir state of the junction is realized. 
The normalization constant is Jo = aDQ 3, /Aed'y 2 kf l . Other 
parameters: DQ 2 — A and the rotation of magnetization in 
the domain walls corresponds to 2Qaq — n. 



small. Therefore, close to the SF interfaces and to the 
boundaries between domains and domain walls (at dis- 
tances of the order of ^) the expressions are not appli- 
cable. 

The main qualitative result of expressions (frjTj) , (p2")) . 
(l70|) . and (|7T|) is that the x component of the current 
density has a form j x = j cx sin(p with negative j cx , i.e., 
the 7T state of the junction is realized. 

The j x {x) function is even while j y (x) is odd with re- 
spect to x = d (the center of the F layer) at both chirali- 
ties. The total current in the y direction is zero. Within 
our accuracy Vj = and V x j = (we must neglect fc w 
and Q in comparison with k n ). 

The LRTC is generated at the boundaries between the 
domains and the domain walls. As a result, the maxi- 
mal supercurrent in the x direction is carried along these 
lines. 

Integrating the x component of the supercurrent 
over y, we find the critical current per period of the 
structured 

r° +QQ / x , aQ 2 nT x ^ /| /fmN 

J c+ = / J CX+ (y)dy = - /I 2 ,4 ElF^' 72 

J-a a -a Q 4e<*tffcJ V K + 

Jc-= j cx -(y)dy = - 2u i l^lr' ( 73 ) 

in the cases of positive and negative chirality, respec- 
tively. The total current across a junction of large area 
is proportional to the number of domain walls. 

The results for the critical current are illustrated in 
Fig. [31 where we plot | J c + | and | J c _ | as functions of the 
domain half-width ap (while Q and oq are linked by the 




FIG. 4: Multidomain F layer of thickness d in contact with a 
bulk superconductor. Depending on the relative orientation 
of rotating magnetizations in neighboring domain walls, we 
distinguish the cases of positive and negative chirality (Q has 
the same or opposite sign in the neighboring domain walls, 
respectively) . The proportion between the widths of domains 
and domain walls is chosen only for drawing purposes. 

condition 2Qciq = n meaning rotation of the magneti- 
zation by the angle n in the domain wall). In the limit 
ao — > the behavior of the supercurrent in the cases of 
positive and negative chiralities is drastically different. 
In the case of positive chirality, the LRTC disappears in 
this limit and the supercurrent vanishes ^ while in the 
case of negative chirality this is the most inhomogeneous 
limit and the supercurrent is maximal. In the opposite 
limit of large ao 3> £t, the ao dependence vanishes in 
both cases and the results coincide since the hyperbolic 
functions of k u ao turn to unity. 

The appearance of the it junction in SFS junctions 
is well understood in the case when it is due to the 
short-range component i The key ingredient is the oscil- 
lating behavior of this component: the additional phase 
7T across the junction is provided by changing the sign. 
At the same time, the LRTC does not change its sign, 
therefore the ir junction due to the LRTC seems counter- 
intuitive. Where does the additional phase come from? 
Note that the LRTC in our case is purely imaginary [see, 
e.g., Eq. ([55]) ]. This means that there is a tt/2 phase rota- 
tion at the SF interfaces, and the two interfaces provide 
the 7T shift. The mechanism of the ir junction due to tt/2 
interface shifts is similar to Ref. 

Another type of SF structures, sensitive to the chirality 
of the vector M, was considered in Refs.[24]and[25|. It was 
shown that the sign of the critical Josephson current in a 
multilayered SF structure depends on chirality. Similarly 
to the present paper, the tt junction was found in the 
situation when only the LRTC is essential. 



IV. DOS IN SF BILAYER 

Our previous paper, Ref. [H, was devoted to studying 
the density of states (DOS) at the free surface of the 
F layer in the system shown in Fig. 0] (with F layer of 
thickness d). These results are immediately reproduced 
from the Green function calculated in the center of the 
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0.05 




FIG. 5: Addition to Fig. 2 from Ref. [23: Correction 5v(y) 
(due to the proximity effect) to the DOS at the free surface 
of the F layer in the case of positive chirality. The curves are 
plotted at several energies e. The width of the domains is 
«o = 5/Q, while the rotation of magnetization in the domain 
walls corresponds to Qaq — it. 



F layer in the SFS junction of Fig. [Tfb) (with F layer of 
thickness 2d) at zero phase difference, <p = 0. Making 
analytical continuation from the Matsubara frequency u> 
to the real energy e, we obtain the correction to the DOS 
due to the proximity effect as 



5v{e) 



(74) 



(we consider the region in space where only the LRTC is 
essential). 

We want to return to this question in view of the recent 
paper Ref. [U where it was demonstrated that general 
analytical properties of the Green function imply that if 
the superconductivity has odd frequency symmetry, then 
5u(0)> 0. The expressions for the DOS from our paper 
Ref. [28] testify that this general statement is satisfied in 
our case, however, this fact was not illustrated in the 
figures. Figures [5] and [6] supplement the figures from 
Ref. [H and demonstrate that 5v(0) > 0. 

Figure [5] is plotted for the same parameters as Fig. 2 
in Ref. |28| and shows the spatial dependence of the DOS 
inside a domain (y is counted from the center of the do- 
main) at several energies. Indeed, the DOS at low ener- 
gies becomes positive everywhere. 

Figure [B] illustrates the 8v{e) dependence at several 
points y. In accordance with Ref. l4ll . the zero-energy 
correction to the DOS is positive forming the zero-energy 
peak. 



V. DISCUSSION 

In our study we neglected orbital effects of the mag- 
netic field in the ferromagnet. Below we demonstrate 
that this is justified for thin F films. We also estimate 
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FIG. 6: Correction 8u(e) to the DOS at the free surface of 
the F layer in the case of positive chirality. The curves are 
plotted at several points y. Other parameters are the same 
as in Fig. [5] 



the influence of the spin-orbit interaction on the LRTC. 
Finally, we comment upon the role of the boundary con- 
ditions at the SF interfaces in our problem. 

The magnetization M in the F layer leads to the ap- 
pearance of the vector potential A which can be esti- 
mated as A ~ AirMd. The vector potential enters Eq. 
([5]) as an additional term (2wA/(j)o) 2 , where 4>o = hc/2e 
is the magnetic flux quantum. This term will restrict the 
penetration length of the LRTC if it is larger than the 
term 2*£ = 2\u\/D ~ 2nT/D. In the opposite limit 



8n 2 Md\ 2ttT 

< 



00 



D 



(75) 



one can neglect the orbital effects (the effect of Meissner 
currents on the LRTC). Taking M ~ 50 G, T ~ 5 K, 
and D ~ 10 cm 2 /s, we obtain d < 300 nm. Therefore 
the orbital effects can be neglected in the case of thin F 
layers. One can also show that under these conditions 
the Meissner currents induced in the superconductors by 
the stray magnetic fields are much smaller than the de- 
pairing currents. Therefore one can neglect the influence 
of the magnetic field of the ferromagnet on the amplitude 
of the order parameter in the superconductors. There is 
one more effect of the domain structure in the ferromag- 
net. This effect occurs in Josephson junctions with lat- 
eral dimensions larger than the Josephson length Aj and 
is related to a spatial variation of the phase difference. 
Due to this a modulation of the total critical Josephson 
current arises. This effect was studied theoretically in 
Ref. S 

If the spin-orbit interaction is present in the F layer, 
it leads to a decrease of the LRTC penetration length. 3 
In this case the wave vector k 2 should be replaced by 

+ A/t so D, where t so is the spin-orbit relaxation time. 
Thus, the LRT C penetration length cannot be larger 
than \J t so D/8. 
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Note also that we use boundary conditions (fT0|) which 
follow from the Kupriyanov-Lukichev ones and describe 
spin-conserving interfaces with potential barriers. These 
boundary conditions themselves do not lead to the ap- 
pearance of the LRTC in the uniform ferromagnet. In 
some papers spin-active interfaces were studied i 42 ' 43 The 
scattering at such interfaces can generate the LRTC even 
in the uniform ferromagnet (see Ref. 44 and references 
therein). The boundary conditions for the spin-active 
interfaces were written in terms of a phenomenological 
scattering matrix. The microscopic calculation of the 
scattering matrix is lacking and the properties of spin- 
active SF interfaces in experiment are not known at 
present. 

In a recent paper Ref. the LRTC was analyzed 
numerically in a model of a SFS structure with a half- 
metallic ferromagnet and spin-active interfaces. The ex- 
change field near the SF interfaces was assumed to be 
inclined with respect to the exchange field in the bulk 
of the ferromagnetic layer. To some extent this model is 
similar to a model of SF structure with a domain wall 
at the SF interface analyzed in Ref. [ill. Therefore, the 
spin-active interfaces are additional sources of the LRTC 
in the ferromagnet. In order to single out the effect of 
the Neel domain walls on the LRTC, we did not consider 
the spin-active interfaces. 

VI. CONCLUSIONS 

We have studied a Josephson junction between two 
superconductors through a multidomain ferromagnet (F) 
with an in-plane magnetization, assuming that the neigh- 
boring domains are separated by the Neel domain walls. 
Due to an inhomogeneous magnetization, the long-range 
triplet superconducting component (LRTC) arises in the 
system. Arising at the domain walls, this component 
spreads into domains over a long distance of the order 
£t = y/D/2TrT, which is much greater than the usual 
short length £/, = y/D/h describing the decaying of su- 
perconducting correlations in a ferromagnet with a ho- 
mogeneous magnetization. 

We have calculated the Josephson current due to this 



component in the case when the short-range components 
exponentially decay over the thickness of the F layer and 
can be neglected. Focusing on the limit when the F layer 
is thin from the viewpoint of the long-range component 
we find that the junction is in the tt state. The LRTC 
does not oscillate inside the F layer, while the additional 
7r phase of the condensate wave function is due to ir/2 
shifts at the SF interfaces. This interpretation suggests 
that the junction must be in the tt state due to the LRTC 
regardless of the F layer's thickness. When the F layer is 
not thin (from the viewpoint of the LRTC), analytical ex- 
pressions for the supercurrent become cumbersome, how- 
ever, numerical calculations indicate that the junction is 
indeed in the tt state. 

The absolute value of the Josephson current density is 
maximal at the boundaries between domains and domain 
walls, see Fig. [2j The current mainly flows along these 
lines. 

We have considered two possible chiralities of the do- 
main structure (positive and negative), which are deter- 
mined by the relative orientation of rotating magneti- 
zations in the neighboring domain walls. The absolute 
value of the Josephson current due to the LRTC is larger 
in the case of the negative chirality, because this case 
corresponds to a more inhomogeneous magnetization. 

Analyzing a correction to the density of states due to 
the LRTC, we find that at zero energy (i.e., at the Fermi 
level) the correction is positive. This fact is in accordance 
with the general statement from Ref. Hi], based on the 
odd frequency dependence of the Green function. 
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